Fluctuation formula for nonreversible dynamics in the thermostated Lorentz gas 
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We investigate numerically the validity of the Gallavotti-Cohen fluctuation formula in the two 
and three dimensional periodic Lorentz gas subjected to constant electric and magnetic fields and 
thermostated by the Gaussian isokinetic thermostat. The magnetic field breaks the time reversal 
symmetry, and by choosing its orientation with respect to the lattice one can have either a generalized 
reversing symmetry or no reversibility at all. Our results indicate that the scaling property described 
by the fluctuation formula may be approximately valid for large fluctuations even in the absence of 
reversibility. 
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PACS numbers: 05.45-a, 05.70.Ln 

The Lorentz gas (LG) thermostated by a Gaussian 
isokinetic (GIK) thermostat is one of the most popular 
models in the study of the relationship between transport 
properties and chaotic behaviour in nonlinear dynamical 
systems. Since the microscopic dynamics of the LG is 
chaotic, and on a sufficiently long time scales it possesses 
a well defined macroscopic transport coefficient, it can be 
used to study the connection of microscopic chaos and 
macroscopic nonequilibrium behaviour. 

The so-called fluctuation formula (FF) has first been 
observed numerically in a system of thermostated fluid 
particles undergoing shear flow 0. In that model, trajec- 
tory segments violating the second law of thermodynam- 
ics were found with probabilities exponentially smaller 
than those of trajectory segments associated with normal 
thermodynamical behaviour. More precisely, let ^r(0 de- 
note the entropy production rate ^ averaged over a time 
interval of length r centered around time t: ^rit) — 

t/2 

- J (,{t + t') dt' , and let us consider it as a probabilistic 

variable. Then its statistical properties in a steady state 
can be characterized by a probability density "Brix). The 
fluctuation formula states P| that this probability density 
has the following property: 



lin, iln^#. 



(1) 



One of the interesting features of the FF is that it seems 
to be valid in systems far from equilibrium, not just for 
vanishing external fields. 

After discovering the formula numerically, analytical 
results were obtained about its validity in deterministic 
systems like transitive Anosov systems ^ and special 
reversible maps 0, Q- In the proofs of these theorems, 
the time reversibility of the system plays a key role |l2j| . 
Nevertheless, proving fluctuation theorems under more 
general conditions seems to be exceedingly difficult. In 
this context, even relatively simple systems like the LG, 



with or without magnetic field, seem to be out of reach 
for the existing analytical techniques. It is also unclear 
how the FF should look like in systems with nonreversible 
dynamics . Consequently, reliable numerical results for 
such models may provide valuable hints in the search for 
more sophisticated theoretical approaches. 

The field driven Lorentz gas consists of a charged par- 
ticle subjected to an electric field moving in the lattice 
of elastic scatterers. For the sake of simplicity, we take 
a square or cubic lattice in our study, depending on the 
dimensionality of the system. Due to the applied elec- 
tric field, one must use a thermostating mechanism to 
achieve a steady state in the system. Such a tool is the 
Gaussian isokinetic thermostat which preserves the ki- 
netic energy of the particle; for a review see e.g. and 
further references therein. We will also apply a constant 
external magentic field to control the reversibility of the 
dynamics. 

Throughout our work we use dimensionless variables. 
We choose the units of mass and electric charge to be 
equal to the mass and electric charge of the particle, so 
we have m = g = 1 in our model. The unit of distance is 
taken to be equal to the radius of scatterers {R= 1), and 
the unit of time is chosen to normalize the magnitude of 
particle velocity to unity. Let q — {qi, . . . ,qn) denote 
the position and p = (pi, . . . ,pn) the momentum of the 
particle in the n-dimcnsional space (n = 2 or 3). Due 
to the normalization, |p| = 1. The phase space variable 
of the system is F = (q, p); it is transformed abruptly 
at every elastic collision and evolved smoothly by the 
differential equation 



q = 
p = 



E + p X B ap 



(2) 
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between them. Here a is called the thermostat variable, 
while E and B are constant vectors playing the roles 
of the external electric and magnetic fields, respectively. 
The GIK thermostat corresponds to the choice a = Ep 
in Eq. For rt = 2, B is thought to be perpendicular 
to the plane of motion given by the directions of E and 
p. We note that Eq. ^ is dissipative, but for B = it 
has also time reversal symmetry. 

Dissipation can be measured by the phase space con- 



2 



traction rate a. It can be computed by taking the diver- 
gence of the right-hand side of Eq. Q : 

cr -divT = - (n - l)a, (3) 

and it can be shown (see e.g. 0) that in our case 

a{t) = m- (4) 

We note that the vahdity of this identity does depend on 
the choosen model and cannot be treated as a general 
property UMM 

The notion of reversibility [ij, an extension of time 
reversal symmetry, can be formulated in terms of the 
phase space flow defined by r{t) = ^'Fq. We say 
that the flow is reversible, if there exists a map G which 
is an involution (i.e is the identity) and bracketing 
the flow by G reverses the direction of time: 



G**G = * 



(5) 



Time reversal symmetry is a special case of reversibility 
with a particular choice of the involution: Go(q, p) — 

(q,-p)- 

In the LG, reversibility depends on the directions of 
the field vectors relative to each other and the lattice. It 
can be checked easily that our system is time reversible 
if B = 0, and it is not otherwise. In Ref. we have 
also shown that the system is still reversible for B 7^ 
if the plane containing E and B is a symmetry plane of 
the lattice. Since the transformation G — MGq (where 
M is a mirroring of q and p with respect to the plane 
containing E and B) satisfies Eq. (O, the smooth flow is 
always reversible. This means that the reversibility of the 
full dynamics including the collisions requires that the 
invariant plane of M be a symmetry plane of the lattice 
[T3I I . In the two dimensional case this is simplified to the 
condition that E has to be contained by the symmetry 
plane of the lattice. 

The goal of our numerical simulations was to measure 
'E,r{x) with a precision which is sufficient to check the 
vahdity of the fluctuation formula. Due to Eq. E.r{x) 
could be measured by periodically computing ar along a 
particle trajectory and making a histogram of these data. 
The disadvantage of this method is that the range of 
possible CTt- values depends on the strength of the electric 
field. Instead we may introduce the quantity 



■Krit) ^ ^ / nEPit + t )dt 



(6) 



where ng denotes the unit vector paralell to E. Since the 
magnitude of p is unity, tt^. always satisfies tt,- G [—1, 1]. 
By making a histogram of the periodically measured val- 
ues of tTt, one gets an approximation of its probability 
density IIt{x). Since the two probabilistic variables sat- 
isfy ar = {n — 1) EtTt, the connection of the probability 
densities is 



1 



(n -l)E \in- l)E 



(7) 
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FIG. 1: The probability density nT(a;) for a 2D configuration, 
where E = (0.5, 0.8) and |B| = 0.2. The distance between the 
centers of the scatterers is d = 2.1; the number of collisions 
is 1.6 X 10*, while the average time between two collisions 
is ~ 0.6. These data are similar throughout all examples 
presented in this paper. 



Then we can rewrite the fiuctuation formula as 
1 



lim 



1 , Tlrix) 

m ■ 



{n-l)ET li^i-x) 



(8) 



At a first glance, 11,- (a;) behaves similarly in all cases: 
as r grows, nT-(a;) becomes more and more concentrated 
around its mean value. This typical shape is shown in 
Fig. It can be noticed that the curve looks like a 
Gaussian, although it is clear that it must be different due 
to the finite range of x 0| . In a separate paper , we 
will deal with the properties of this distribution in more 
details. In order to visualize the fluctuation formula, we 
introduce the quantity 



Drix) 



1 



l^^UAx) 



{n-l)ET Uri-x) 



(9) 



that must exactly be linear with a slope 1 in the t — > cxd 
limit if the fluctuation formula is valid. We will investi- 
gate for different configurations of the LG how well Dr{x) 
approaches this behavior in numerical simulations. Due 
to the fact that we have a finite number of data points 
coming from a numerical trajectory of finite length, our 
conclusions concerning Dr{x) and thus the fluctuation 
formula are, of course, limited to an interval [— A^, A,-] 
with A^ < 1. In practice, if the extremal tTt values in a 
series of A'" measured data were TTmin and iTmax , then we 
identified A^. with min(— 7r„jj„, i^max)- It is easy to check 
that the probability of observing tt^. values outside this 
interval in another series is in the order of l/N. Fig. |21 
shows the r dependence of A,- for different field strengths 
with N fixed. For our simulations, N was chosen to be 
10^, which means that if Dt-{x) is found to be linear on 
[—A,-, A7.] with slope 1, then it can be interpreted as the 
fluctuation formula is valid for fluctuations with proba- 
bilities larger than 10~^. 

In the rest of the paper, we present our numerical re- 
sults for the GIK thermostatted LG both in two and three 
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FIG. 2: The dependence of A,- on r for diflerent field 
strengths in a 2D configuration with B = 0. The direction of 
E is parallel with (5, 8) but its magnitude varies. The curves 
appear to be linear in the dominant region on the log- log plot, 
suggesting a power law dependence on r. This behavior seems 
to be valid for other configurations as well, no matter they are 
reversible or not. 
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FIG. 3: Dr(x) for a time reversible configuration (B = 0) 
in 2D with E = (0.5, 0.8). The inset shows that for higher r 
values, Dt{x) a; on [— Ar, A^]. The inset has the same axes 
as the figure. 
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FIG. 4: Dr(x) for a time reversible configuration in 3D, with 
E = (0.05, 0.1, 0.15) and B = 0. The inset shows that Dr{x) 
converge to a; as r gets larger. The axes of the inset are the 
same as in the figure. 
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FIG. 5: D-r{x) for a nonreversible configuration in 2D, with 
|B| = 0.2 and E = (0.5, 0.8). It seems that for lower r values 
D-rix) has a breakoff from the linear curve around x ~ ±0.3, 
but the inset shows that for higher r values, Dt{x) behaves 
quite similarly to the reversible case of Fig. El We note that 
this is the same configuration as the one used for Fig. The 
inset has the same axes as the figure. 



dimensions, with various values of the external fields. We 
focus on the question whether nonreversible dynamics 
leads to different scaling in the fluctuations than the one 
found in reversible systems. As a general rule, we have 
not found any difference between time reversal symmet- 
ric cases (i.e., with B = 0) and reversible ones. Indeed, 
time reversal symmetry can be replaced in the known 
fluctuation theorems by general reversibility without af- 
fecting their validity, since the proofs do not make use 
of the special form of the involution Go- We note that 
we have tested the different dynamical cases with several 
choices for the field stregths and could not find signifi- 
cant deviations in the observed behavior as long as we 
stayed within the ergodic region of the parameter space. 

In Figs. 121 and ^ we plot Dr {x) for reversible dynam- 
ics in two and three dimensions (2D and 3D), respec- 
tively. The fiuctuation formula appears to be valid in 



both cases; the convergence to the linear limit, however, 
seems to be different in them. For 2D, the Dt {x) curve 
has deviations, decreasing in size with t increasing, from 
the linear shape, while for 3D, Dr {x) exhibits strongly 
linear behavior with slopes approaching 1 as t increases. 
It is worth noting that the latter convergence can also be 
observed in the 2D random LG • 

Our results for the nonreversible versions are shown in 
Figs. O and O The most striking difference compared 
to the reversible cases is the fact that there seems to be 
a cubic term present in Dr {x) that does not disappear 
for larger r values. This term leads to a breakoff from 
the diagonal line for \x\ > Xc ~ 0.3, which means that 
there can be deviations from the fluctuation formula for 
large fluctuations. The slope of the linear part, how- 
ever, is still 1 in the large r limit, so the FF can be a 
good approximation for small to moderate size fluctu- 
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FIG. 6: Dt(x) for a nonreversible configuration in 3D, with 
E = (0.05,0.1,0.15) and B = (0.16, -0.06,0.04). The axes of 
the inset are the same as in the figure. 



more and more valid in a probabilistic sense as r — > oo, 
since the larger fluctuations become less and less likely 
in that limit. 

We may conclude that the FF appears to be valid in the 
GIK thermostated LG with reversible dynamics, both in 
two and three dimensions. For nonreversible dynamics, 
we have found indications that the FF may still describe 
the scaling properties of fluctuations in a moderate size 
regime, although for large fluctuations there are clear de- 
viations from it due to higher order terms in Dt {x) . The 
fact that the slope of the linear part in the scaling behav- 
ior is the same in reversible and nonreversible cases sug- 
gests a kind of robustness for the FF in the thermostated 
LG. It would be interesting to see if this remains valid in 
other nonequlibrium systems with nonreversible dynam- 
ics. 



ations. The fact that the region of validity of the FF 
does not shrink considerably for larger t values suggests 
that the coefficient of the cubic term in Dt-{x) may have 
only weak dependence on r. This also means that as 
the distribution T\r{x) is concentrating around its mean 
value for increasing r values, the total statistical weight 
of the large fluctuations that are not covered by the lin- 
ear regime is decreasing. In other words, the FF becomes 
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